FAST CARS: Engineering a Laser Spectroscopic Technique for 
Rapid Identification of Bacterial Spores 

M. O. Scully,* G. W. Kattawar/'^ R. P. Lucht/'^ T. Opatrny/'^ H. Pilloff/ 
A. Rebane/ A. V. Sokolov/'^ and M. S. Zubairy^-^'^ 
Q I Institute for Quantum Studies^, Dept. of Physics^ , Dept. of Electrical Engineering^ , Dept. 

^ ■ of Mechanical Engineering'^, Texas A&M University, College Station, Texas 77843 

^1 ^ Max- Planck- Institut fiir Quantenoptik, D-85748 Garching, Germany 



(N 



^Dept. of Theoretical Physics, Palacky University, Olomouc, Gzech Republic 
Dept. of Physics, Montana State University, Bozeman, Montana 59715, USA 



\ ^Dept. of Electronics, Quaid-i-Azam University, Islamabad, Pakistan 

'• (February 2, 2008) 

q ; 

c/5 . Airborne contaminants, e.g., bacterial spores, are usually analyzed by time consuming 

' microscopic, chemical and biological assays. Current research into real time laser spectro- 

^ 1 scopic detectors of such contaminants is based on e.g. resonance fluorescence. The present 

Oh! approach derives from recent experiments in which atoms and molecules are prepared by one 

^ . (or more) coherent laser(s) and probed by another set of lasers. These studies have yielded 

^> ' such counterintuitive results as lasers which operate without inversion, ultra-slow light with 
t — ■ 

O i group velocities of order 10 meters/sec, and generation of ultra-short pulses of light via 

^ ! phased molecular states. The preceding examples are based on inducing a phase coherent 

^ . state of matter in the ensemble of simple molecules being studied. The connection with 

O . previous studies based on "Coherent Anti-Stokes Raman Spectroscopy" (CARS) is to be 

Q ' noted. However generating and utilizing maximally coherent oscillation in macromolecules 

^ I having an enormous number of degrees of freedom is much more challenging. In particular, 

^ . the short dephasing times and rapid internal conversion rates are major obstacles. However, 

^ ' adiabatic fast passage techniques and the ability to generate combs of phase coherent fem- 

' tosecond pulses, provide new tools for the generation and utilization of maximal quantum 

^ I coherence in large molecules and biopolymers. This extension of the CARS technique is 
called FAST CARS (Femtosecond Adaptive Spectroscopic Techniques for Coherent Anti- 
Stokes Raman Spectroscopy), and the present paper proposes and analyses ways in which 
it could be used to rapidly identify pre-selected molecules in real time. 



My friend and mentor Vicky Weisskopf used to say "The best way into a new problem is to bother people." 
This is faster than searching the literature and more fun. I would like to thank my colleagues for allowing 
me to be a bother and especially my coauthors who have suffered the most! This paper is dedicated to 
Prof. Viktor von Weisskopf: premier physicist and scientist-soldier who stood by his adopted country in 
her hour of need. 
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I. INTRODUCTION 

There is an urgent need for the rapid assay of chemical and biological unknowns, such 
as bioaerosols. Substantial progress toward this goal has been made over the past decade. 
Techniques such as fluorescence spectroscopy [|l], |[, and UV resonant Raman spectroscopy 
0, ^, H, H have been successfully applied to the identification of biopolymers, bacteria, and 
bioaerosols. 

At present field devices are being engineered [|T| which will involve an optical preselection 
stage based on, e.g., fluorescence radiation as in Fig. |^. If the fluorescence measurement 
does not give the proper signature then that particle is ignored. Most of the time the 
particle will be an uninteresting dust particle; however, when a signature match is recorded, 
then the particle is selected for special biological assay, see Fig. Pd. The relatively simple 
fluorescence stage can very quickly sort out some of the uninteresting particles while the 
more time consuming bio-tests will only be used for the "suspects" . 

The good news about the resonance fluorescence technique is that it is fast and simple. 
The bad news is that while it can tell the difference between dust and bacterial spores, it 
can not differentiate between spores and many other organic bioaerosols, see Fig. |l|c. 

However, in spite of the encouraging success of the above mentioned studies, there is 
still interest in other approaches to, and tools for, the rapid identification of chemical and 
biological substances. To quote from a recent study [0: 

"Current [fluorescence based] prototypes are a large improvement over earlier 
stand-off systems, but they cannot yet consistently identify specific organisms 
because of the similarity of their emission spectra. Advanced signal processing 
techniques may improve identification." 

Resonant Raman spectra hold promise for being spore specific as indicated in Fig. §b. 
This is the good news, the bad news is that the Raman signal is weak and it takes several 
minutes to collect the data of Fig. ^3. Since the through-put in a set-up such as that of 
Fig. |l|b is large, the optical interrogation per particle must be essentially instantaneous. 

The question then is: Can we increase the resonant Raman signal strength and thereby 
reduce the interrogation time per particle? If so, then the technique may also be useful in 
various detection scenarios. 
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The answer to the question of the proceeding paragraph is a quahfied "yes." We can 
enhance the Raman signal by increasing the coherent molecular oscillation amplitude Rq in- 
dicated in He. In essence this means maximizing the quantum coherence between vibrational 
states \b) and |c) of Qa. 

Our point of view derives from research in the fields of laser physics and quantum optics 
which have concentrated on the utilization and maximization of quantum coherence. The 
essence of these studies is the observation that an ensemble of atoms or molecules in a 
coherent superposition of states represents, in a real sense, a new state of matter aptly 
called "phaseonium" [§]. 

In particular, we note that matter in thermodynamic equilibrium has no phase coherence 
between the electrons in the molecules making up the ensemble. This is discussed in detail 
in Section III. When a coherent superposition of quantum states is involved, things are very 




FIG. 1: (a) Ultra-violet (UV) excitation radiation promotes molecules from ground state 
to an excited state manifold. This excited state manifold decays to the ground state 
via non-radiative processes to a lower manifold which then decays via visible or UV flu- 
orescence. It is this fluorescence signature which is detected in resonance fluorescence 
spectroscopy; (b) Figure depicting a scenario in which a UV laser interacts with dust 
particles and bio-spheres of interest. When, for example, a bacterial spore is irradiated, 
fluorescence will be emitted signaling that this particular system is to be further tested. In 
principle, uninteresting particles are deflected one way; but when fluorescence takes place, 
the particles are deflected in another direction and these particles are then subjected to 
further biological tests; (c) The shaded area displays the signal range for the fluorescence 
spectrum of a number of biological samples. Bacillus subtilis, Bacillus thuringiensis, Es- 
cherichia coli, and Staphylococcus aureus. It is not possible to distinguish between the 
different samples based on such a measurement (see Q for more details). 
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FIG. 2: (a) Resonant Raman scattering in which radiation ui excites the atom 
from |c) to I a) and the Stokes radiation is emitted taking the molecule from \a) 
to 1 6). The frequency of the excitation radiation is i^i, and the frequency of the 
Stokes radiation is 1^2', (b) Detail of UV resonance Raman spectra of spores of 
Bacillus megaterium (1), Bacillus cereus (2), and of Calcium Dipicolinate (3), all 
excited at 242 nm; adapted from W.H. Nelson and J.F. Sperry, Modern Techniques 
in Rapid Microorganism Analysis, edited by W.H. Nelson (VCH Publishers, N.Y. 
1991). P (see also Fig. (c) Figure illustrating a more physical picture of Raman 
scattering in which a single diatomic molecule, consisting of a heavy nucleus e.g. 
Uranium 235 and a light atom e.g. hydrogen, scattering incident laser radiation 
at frequency vi. The vibrational degrees of freedom associated with the diatomic 
molecule are depicted here as occurring with amplitude Rq oscillating at frequency lo. 
The scattered radiation from this vibrating molecule is at frequency 1^2 = ^^i — "-^ for 
the Stokes radiation. This classical picture of the vibrating dipole is to be understood 
as an amplitude times a sinusoidal oscillation at the frequency uj as indicated in the 
equation beneath the figure. Rq is a quantum mechanically calculated oscillation 
amplitude as discussed in detail in Section III of the paper. 



different and based on these observations, many interesting and counterintuitive notions 
are now a laboratory reality. These include lasing without inversion (LWI) 0, electromag- 
netically induced transparency (EIT) light having ultra slow group velocities on the 



order of 10 meters/sec |TT]], and the generation of ultra short pulses of light based on phased 
molecular states [O. 
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Another emerging technology central to the present paper is the exciting progress in the 
area of femtosecond quantum control of molecular dynamics originally suggested by Judson 



and Rabitz [0. This is described and reviewed in the articles by Kosloff et al. [|T^], Warren, 
Rabitz and Dahleh |T^, Gordon and Rice |T^, Zare Rabitz, de-Vivie-Riedle, Motzkus 
and Kompa |T8[, and Brixner, Damrauer and Gerber ||19|. Other related work on quantum 



coherent control includes: The quantum interference approach of Brumer and Shapiro P 



the time-domain (pump-dump) technique proposed by Tannor, Kosloff and Rice the 



stimulated Raman Adiabatic Passage (STIRAP) approach of Bergmann and co-workers |22 
to generate a train of coherent laser pulses. The preceding studies teach us how to produce 
pulses having arbitrary controllable amplitude and frequency time dependence. Indeed the 
ability to sculpt pulses by the femtosecond pulse shaper provides an important new tool 
for all of optics, see the pioneering works by Heritage, Weiner, and Thurston |^|, Weiner, 
Heritage, and Kirschner |2^, Wefers and Nelson p5| and Weiner p6[| . 

An important aspect of the learning algorithm approach is that knowledge of the molec- 
ular potential energy surfaces and matrix elements between surfaces are not needed. Precise 
taxonomic marker frequencies may not be known a priori; however, by using a pulse shaper 
coupled with a feedback system, complex spectra can be revealed. 

Thus, we now have techniques at hand for controlling trains of phase coherent femtosecond 
pulses so as to maximize molecular coherence. This allows us to increase the Raman signal 
while decreasing the undesirable fluorescence background. This has much in common with 
the CARS spectroscopy of Fig. ||, but with essential differences as we now discuss. 

The presently envisioned improvement over ordinary CARS is based on enhancing the 
ground state molecular coherence. However, we note that molecules involving a large number 
of degrees of freedom will quickly dissipate the molecular coherence amongst these degrees 
of freedom. This is a well known difficulty and is addressed in the present work from several 
perspectives. First of all, when working with ultra short pulses, we have the ability to gen- 
erate the coherence on a time scale which is small compared with the molecular relaxation 
time. Furthermore, we are able to tailor the pulse sequence in such a way as to mitigate, and 
overcome key limitations in the application of conventional CARS to trace contaminants. 
The key point is that we are trying to induce maximal ground state coherence, as opposed 
to the usual situation within conventional CARS where the ground state coherence is not a 
maximum as is shown later in this paper. With FAST CARS (Femtosecond Adaptive Spec- 
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troscopic Techniques applied to Coherent Anti-Stokes Raman Spectroscopy) we can prepare 
the coherence between two vibrational states of a molecule with one set of laser pulses; and 
use higher frequency visible or ultra-violet to probe this coherence in a coherent Raman 
configuration. This will allow us to capitalize on the fact that maximally coherent Raman 
spectroscopy is orders of magnitude more sensitive than incoherent Raman spectroscopy. 

Having stated our goals and our approach toward attaining these goals, we emphasize 
that the present paper represents essentially an engineering endeavor. We propose to draw 
heavily on the ongoing work in quantum coherence and quantum control as mentioned 
earlier. 

For example, the careful experiments and analysis of the Wiirzburg group on the gen- 
eration and probing of ground state coherence in porphyrin molecules by femtosecond- 
CARS (fs-CARS) are very germane to our considerations. However, ground state coherence 
is not maximized in these experiments. 



(a) 
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FIG. 3: (a) Ordinary resonant Raman spectroscopy in which a drive laser of amplitude 
£i generates a weak signal field having an amplitude £2- The incident signal consists of 
one pulse at vi and the pulse structure following interaction with the molecular medium 
consists of two pulses at i^i and 1^2- (b) The coherent Raman process associated with 
CARS is depicted in which two fields at frequency i^i and 1^2 are incident with amplitudes 
£1 and £2- The third radiated anti-Stokes signal field at frequency 1^3 is indicated. Hence 
CARS involves 2 fields in and 3 out. (c) FAST CARS configuration in which maximal 
coherent Raman spectroscopy is envisioned. The preparation pulses £1 and £2 prepare 
maximum coherence between states |6) and |c). Next the probe laser £3 interacts with 
this oscillating molecular configuration and the anti-Stokes radiation is generated. Thus 
we have 3 fields in and 4 out when using FAST CARS. 
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In another set of beautiful experiments ||2^ they investigate the selective excitation of 
polymers of diacetylene via fs-CARS. They control the timing, phase and frequency (chirp) 
content of their preparation pulses. In these experiments it was necessary to focus attention 
on the evolution of the excited state molecular dynamics. We hope to avoid this complication 
as is explained later. 

Perhaps closest to our approach is the recent joint work of the Garching Max-Planck 
and Wiirzburg groups |^ . Their paper entitled "Optimal control of ground-state dynamics 
in polymers" is a prime example of a FAST CARS experiment. However they concentrate 
on producing highly excited states of the "vibrational motion of a certain bond". The 
application of their technique to the production of maximum coherence between states \b) 
and |c) of Fig. ^a in a specific vibrational mode of their molecule would be of great interest 
to us and is underway. 

Finally we wish to draw the reader's attention to the useful collection of articles in a 



recent special issue of the "Journal of Raman Spectroscopy" dedicated to fs-CARS |^ 



Likewise the recent work of Silberberg and coworkers in which they show that it is 
possible to excite one of two nearby Raman levels, even when they are well within the broad 
fs pulse spectrum is another excellent example of the power of the FAST CARS technique. 

To summarize: the present work focuses on utilization of a maximally phase coherent 
ensemble of molecules, i.e. molecular phaseonium, to enhance Raman signatures. This 
will be accomplished via the careful tailoring of a coherent pulse designed to prepare the 
molecule with maximal ground state coherence. Such a pulse is a sort of "melody" designed 
to prepare a particular molecule. Once we know this molecular melody, we can use it to set 
that particular molecule in motion and this oscillatory motion is then detected by another 
pulse; this is the FAST CARS protocol depicted in Figs ^ and p!5|b. 

In order to establish the viability and credibility of this program, the material covered 
in the present paper is presented in some detail. It is hoped that scientists who are experts 
in one phase of the subject, e.g., molecular biology but not with subtleties of modern laser 
spectroscopy can read the paper without undue appeal to the literature or complicated 
mathematical developments. On the other hand, some basic facts of life, endosporewise, are 
important. Hence, a short overview of some aspects of Raman spectroscopy as applied to 
macromolecules and especially to biological spores is presented. 

In Section II, the status of Raman spectroscopy applied to biological spores is reviewed. 
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In Section III, we compare various types of Raman spectroscopy with an eye to the recent 
successful apphcations of quantum coherence in laser physics and quantum optics. Section 
IV presents several experimental schemes for applying these considerations to the rapid 
identification of macromolecules, in general, and biological spores, in particular. Finally in 
Section V we propose several scenarios in which FAST CARS could be useful in the rapid 
detection of bacterial spores. Where appropriate, mathematical details are included in 
Appendices and comparison between the various types of Raman spectroscopic techniques 
are discussed with special emphasis on overall sensitivity. As stated earlier, the present 
paper is an engineering science analysis of a promising approach to the problem of bacterial 
spore detection. This is not a review paper. If the reader feels that we have missed or 
misrepresented her research, we would be happy to learn how to use it to improve our 
design and detection strategy. Indeed, we view this paper as providing a point of departure, 
and will be pleased if it provokes discussion and debate. If the reader is not provoked we 
apologize. It is very difficult to annoy everybody in a single paper. 
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II. PICO-REVIEW OF RAMAN SPECTROSCOPY APPLIED TO BACTERIAL 

SPORES 



The bacterial spore is an amazing life form. Spores thousands of years old have been 
found to be viable. One textbook ||3^ reports that "endospores trapped in amber for 25 
million years germinate when placed in nutrient media." 

A key to this incredible longevity is the presence of dipicolinic acid (DPA) and its salt 
calcium dipicolinate in the living core which contains the DNA, RNA, and protein as shown 
in Fig. g 
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FIG. 4: Sketch of spore indicating that the DPA and its salts, e.g., Ca-DPA are contained 
in the core and are in contact with the spore specific DNA ribosomes and cell proteins. 

A major role of the calcium DPA complex seems to be the removal of water, as per the 
following quote |]3^ : "The exact role of these [DNA] chemicals is not yet clear. We know, 
for instance, that heat destroys cells by inactivating proteins and DNA and that this process 
requires a certain amount of water. Since the deposition of calcium dipicolinate in the spore 
removes water ... it will be less vulnerable to heat." 

Hence, one of the major components of bacterial spores is dipicolinic acid (DPA) and its 
ion as depicted in Fig. |^. Calcium dipicolinate can contribute up to 17% of the dry weight of 
the spores. A definitive demonstration of this conjecture was made by comparing the 242 
nm excitation spectra of calcium dipicolinate with spore suspensions of Bacillus megaterium 
and Bacillus cereus. From Fig. |^, it is seen that good matches were noted for the 1017, 
1396, 1446, and 1607 cm^^ peaks of the calcium dipicolinate. 

As has been emphasized by W. Nelson and coworkers @, ^ 0, the presence of DPA 
and its calcium salt gives us a ready made marker for endospores. As has been mentioned 
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FIG. 5: (a) Dipicolinic acid (2, 6-pyridinedicarboxylic acid, C5H3N(COOH)2); (b) The 
Ca^+ DPA complex. 



earlier and as will be further discussed later, this is the key to Raman fingerprinting of the 
spore. 

We note however that fluorescence spectroscopy was one of the first methods used for 
detection of bacterial taxonomic markers and is still used for detection where high specificity 
is not required. This technique is an important addition to the "tool kit" of scientists and 
engineers working in this area. 

1046 
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FIG. 6: This figure (adapted from shows UV resonance Raman spectra of spores of 
Bacillus megaterium (1), spores of Bacillus cereus (2), and calcium dipicolinate (3) in 
three spectral regions. All samples are excited at 242 nm. 

A possible FAST CARS protocol is as follows: First we obtain size and fluorescence 
information. If this is consistent with the presence of a particular bacterial spore we could 
then automatically perform a FAST CARS analysis sensitive to DPA so as to further narrow 
the number of suspects. 
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It is important to note that just as resonant Raman is some 10^ times more sensitive 
than non- resonant, coherent Raman yields a much stronger signal than ordinary incoherent 
Raman spectroscopy. This makes it possible to collect the Raman spectra much more rapidly 
via FAST CARS and this is very important in the ultimate scheme of things. 

To summarize: we will generate quantum coherence in macromolecules by working with 
the now available femtosecond pulse trains in which there exists phase coherence between 
the individual pulses. In this way, one can enhance coherent Raman signatures. The uti- 
lization of "molecular music" to generate maximal phase coherence holds promise for the 
identification and characterization of macro and bio-molecules. 
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III. COMPARISON OF DIFFERENT TYPES OF RAMAN SPECTROSCOPY 



Raman scattering is an inelastic scattering of electromagnetic fields off vibrating 
molecules. The origin of Raman scattering dates back to a theoretical paper in Naturwis- 
senschaften by A. Smekal in 1923 entitled (translated) "The quantum theory of dispersion" 



35|. It was followed by another paper in a 1923 Physical Review (by A. Compton) entitled 



"A quantum theory of the scattering of X-rays by light elements" ||36|. Some historians feel 
that these two papers gave CV. Raman the idea for the experiments that were performed 
with K.S. Krishnan and led to the discovery of the effect in over 60 liquids. Raman and 
Krishnan published their results entitled "A new type of secondary radiation" in Nature 
on March 28, 1928 |3^. It was soon followed by the landmark paper of G. Landsberg and 
L. Mandelstam who found the same effect in quartz and published a paper entitled (trans- 
lated) 'A novel effect of light scattering in crystals" which appeared on July 13, 1928 in 
Naturwissenschaften [^. By the end of 1928 dozens of papers had already been published 
on the "Raman" effect. 

In this section we first recall the quantum mechanical picture of a vibrating molecule. 
We then discuss the principles of different types of Raman spectroscopy. 



A. Molecular vibrations 



Let us consider a simple diatomic molecule for explanation of the principle. The in- 
teratomic oscillation can be visualized via a classical picture of the vibrating molecule as 
in Fig. Quantum mechanically, the situation can be understood as depicted in Fig. 
^. The potential energy of the molecule depends on the interatomic distance z and has a 
well pronounced minimum. The Hamiltonian of the vibrating molecule has a set of dis- 
crete eigenstates; in Fig. ^ we show just the ground state |c) and the first excited state \h). 
Whereas in each of these states the mean displacement from the equilibrium position is zero, 
a quantum superposition of these states has generally a nonzero mean displacement R{t) 
which varies with time. Assuming that in time t = the molecule is in a superposition state 
|\l/(0)) = b\b) + c|c), then in time t > the state is |\I/(t)) = b\b) + cexp{—iu!t)\c) , where 
the frequency u = {Ec — Eb)/h is the difference of the energies of the eigenstates |c) and \b) 
divided by the Planck constant h. The mean displacement R{t) is then R{t) = (t)\R\'^ (t)) 
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is then 

R{t) = Rq exp{ujt — (Pq) + c.c, (1) 

where 

Ro = \bc*\{b\R\c) (2) 

is the displacement amphtude, ipo is the initial phase (determined by the phases of the 
coefficients b and c), and R is the displacement operator. It can be seen from Eq. 
that one can reach the maximum amplitude of the mean displacement if the superposition 
coefficients b and c are of the same mae: nitude, i.e., \b\ = \c\ = 1/V2 so that = 
{b\R\c)/2. Thus the product be* is of special importance for determining the vibrational 
amplitude. This coherent superposition of states is generally described by the off-diagonal 
density matrix element pbc which for the present simple case is given by = be*. Without 
going into detail we simply state that the density matrix element pf,c is a complex number 
(0 < \pbc\ < 1/2) characterizing the quantum state of the molecule and determining the 
amplitude of the mean displacement. For some quantum states quantum coherence is not 
present (e.g., energy eigenstates, thermal states, etc.), whereas for some states it can reach 
the maximum magnitude (i.e., for |\E') = 2~^/^(|6) + |c))). 

We emphasize that when all molecules are in the same superposition state, the response 
of the sample to an optical signal is very different from the thermal state. Preparation and 
optical probing of molecular vibrations is the essence of Raman spectroscopy. 

B. Classical description of the Raman scattering 

The simplest classical description of the Raman effect assumes that the polarizability a of 
the molecule is dependent on the relative positions of the atomic nuclei. The polarizability 
is the proportionality factor between the external electric field E and the molecular dipole 
moment P, P = aE, and for a vibrating molecule it is a time-dependent quantity. In the 
linear approximation, and assuming just the scalar case, the polarizability can be written as 
a{Q) = a(0) + a'Q, where Q is the generalized coordinate of the vibrating molecule, a(0) 
is the polarizability of the equilibrium state, and a' = da/dQ with Q = 0. If the molecule 
vibrates with frequency Ubc, the coordinate Q changes as Q = Qo cos^Uabt)- The electric field 
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Coherent: pUt) = b(t)c*(t) Thermal: Pbc= ^ 



FIG. 7: Molecular coherence: atomic distance of a molecule in a coherent superposition 
of the vibrational states \b) and |c) oscillates periodically with time. On the other hand, 
a molecule in a thermal equilibrium state is stationary and its coherence g^c vanishes. 

irradiating the molecule oscillates as E — Eq cos(i/ii) . Thus, one finds that the molecular 
dipole oscillates with several frequencies: with the frequency of the incoming radiation ui 
(leading to the Rayleigh scattering), and with the shifted frequencies ui ± uj^c (leading to 
the Stokes and anti-Stokes Raman frequencies). 

Even though this model is able to predict the correct frequencies of the scattered fight, 
it cannot tell us anything about the intensities of different field components in spontaneous 
scattering processes. To get more information about the scattering process, one needs a 
quantum mechanical model of the molecule. Let us now study the main features of the 
various Raman scattering processes. 

C. Stokes vs. anti-Stokes scattering 

Raman scattering is an optical phenomenon in which there is a change of frequency of 
the incident light. Light with frequency ui scatters inelastically off the vibrating molecules 
such that the scattered field has frequency 1/2 — J^i^ i^bc, where uj^c is the frequency of the 
molecular vibrations. The field with down-shifted frequency 1/2 — J^i — oJbc is called Stokes 
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field and its generation corresponds to the process depicted in Fig. |a, whereas the frequency 
up-shifted radiation is called the anti-Stokes field and corresponds to the process in Fig. |b. 



(a) (b) 
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FIG. 8: Stokes (a) and anti-Stokes (b) Raman scattering. Pump field 1 interacts with 
a vibrating molecule to produce the scattered field 2 which has either lower frequency 
(Stokes scattering) or higher frequency (anti-Stokes scattering). 



D. Spontaneous vs. stimulated Raman scattering 

There are two basic Raman processes: the so-called spontaneous and stimulated Raman 
scattering. Spontaneous scattering occurs if a single laser beam with intensity below a 
certain threshold illuminates the sample. In condensed matter, in propagating through 1 
cm of the scattering medium, only approximately 10~^ of the incident radiation is typically 
scattered into the Stokes field (see, e.g., ||3^). Stimulated scattering which occurs with a 
very intense illuminating beam is a much stronger process in which several percent of the 
incident laser beam can be converted into the other frequencies. From the quantum-optical 
point of view, the Raman scattering can be described by means of photon numbers occupied 
in different modes. The rate of photon number increase in the Stokes mode can be written 
as ris = ^niins + 1), where ni is the number of photons in the incident laser mode and ns 
is the number of photons in the Stokes mode. Here ?7 is a proportionality constant. 

For spontaneous Stokes scattering 71,5 <^ 1, and the intensity of the scattered field is 
roughly proportional to the length traveled by the incident field in the medium. On the 
other hand, for 7^5 > 1, the stimulated process becomes dominant and the scattered field 
intensity can increase exponentially with the medium length. 
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TABLE I: Comparison of different Raman spectroscopic techniques as derived in Appendix The 
density matrix element pbc governs the amphtude of coherent vibration, whereas the element pab 
is proportional to the electronic polarization responsible for emission of radiation. r2i^2 are the 
Rabi frequencies, A is the detuning of the electronic transition, ^yab, Jac are the decay rates of 
the optical transitions, jbc is the decoherence rate of the vibrational states, and 71 is the decay 
rate from level b to c. The approximated values (shown in the lower right corner) were obtained 



for Jab ~ la 



Ibc 



12^-1 
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lO^s ^, A w lO-'^^s ^, and r2i_2 ~ lO-'^^s ^ for weak driving 



and r2i^2 ~ 10 s for strong driving. Note that Q. ~ 10 s corresponds to a 10 ns pulse with 
0.1 mJ energy focused on a square millimeter spot if the electronic transition dipole moment is 
p « 10-i9CxlO~^°m [see Eq. (|A1^)]. 
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E. Resonant vs. non-resonant Raman processes 

The resonant Raman process (appearing when the frequency of the incident radiation 
coincides with one of the electronic transitions) is much richer than the nonresonant, and 
we now turn to a discussion of the resonant problem. 

Resonant Raman radiation is governed by the oscillating dipole between states \a) and 
\b) (Stokes) and/or \a) and |c) (anti-Stokes) in the notation of Fig. |^ and Table |. In the 
Stokes case, the steady state coherent oscillating dipole P(t), divided by the dipole matrix 
element pab = e{a\r\b), is the important quantity. That is pab{t) = P{t)/pab, as given by 
Eq. (^, is 

Qab = -i [^2{na - Ub) - VtiQcb] / Ylab " «(t^afe - 1^2)] (3) 

where the Raman coherence is pbc as discussed earlier, e.g. Fig. 0. In Eq. (|]), ojab is 
the transition frequency between the electronic states a and 6, z/2 is the frequency of the 
generated field, and the other quantities are defined in the caption of Table |. 

The main advantage of resonant Raman scattering is that the signal is very strong — up 
to a million times stronger compared to the signal of nonresonant scattering [^]. It is also 
very useful that only those Raman lines corresponding to very few vibrational modes associ- 
ated with strongly absorbing locations of a molecule show this huge intensity enhancement. 
On the other hand, the resonance Raman spectra may be contaminated with fluorescence. 
However, this problem can be avoided by using UV light so that most of the fluorescence 
appears at much longer wavelengths than the Raman scattered light and is easily filtered 
out. 

F. Coherent vs. incoherent Raman scattering from many molecules 

An important distinction between different Raman scattering schemes is based on the 
phase relation of the field scattered off different molecules. In the incoherent case the spon- 
taneous contributions of individual molecules sum up with random phases. The magnitude 
of the emitted electric field then scales as y/N (see Fig. |^). On the other hand, if all 
the molecules are prepared in the same coherent superposition of their states \h) and |c), 
their contributions to the emitted field have the same phases and the magnitude of their 
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FIG. 9: Incoherent (a) vs. coherent (6) addition of the radiation from N molecules. Each 
arrow represents the contribution of one molecule to the resulting electric field vector. 



sum is proportional to (see Fig. ^b). Thus, the intensity radiated from the sample is 
proportional to in the incoherent case and to A^^ in the coherent case. 

There is another important aspect of the coherent vs. incoherent resonant scattering 
process, namely the rate of emission from the atom ensemble. In the far off resonance, 
Smekel-Raman limit, the emission and absorption are simultaneous since the transitions to 
the excited state(s) are virtual. This is not true for the resonant Raman processes. In that 



case the molecule is excited to the |a) state of Fig. |T0| where, in the case of weak Stokes 
field it can live for many nanoseconds. 

However, the cooperative emission rate from A^ molecules can be much faster than that 
of a single molecule. As was shown by Dicke, and derived in Appendix 0, if the spontaneous 
emission lifetime of A^ incoherent dipoles is given by r; then the cooperative spontaneous 
emission lifetime of A^ coherently prepared dipoles can be as short as r/A^. 

This superradiant "speed up" of the radiation process can have important consequences 
for the present problem. We recall that the experiments of Nelson and coworkers p, |, |^, |^ 
are carried out at 242 nm so as to "ride above" the fluorescence noise, see Fig. p^b. But, 
if we can enhance the spontaneous emission rate so that the transition rate from the a2 
manifold to the ground state g is faster than from the internal (non radiative) rate from 
02 — i> ai, it would be possible to mitigate fluorescence noise. Then it would be possible to 
carry out resonant Raman studies with visible or near UV lasers instead of using the 242 
nm wavelength. 
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IV. FAST CARS 



A. Generation of atomic coherence 



The purpose of this section is to demonstrate the utihty of pulse shaping as a mechanism 
for generating maximal coherence. The Raman signal is optimized at the condition of 
maximal molecular coherence. When in this state, each of the molecules oscillates at a 
maximal amplitude, and all molecules in an ensemble oscillate in unison. Here we discuss 
several methods for the preparation of maximal coherence state. 



1. Adiabatic Rapid Passage via Chirped Pulses 



A particularly simple and robust approach to the generation of the maximal coherence 
is to use a detuning 6uj which is largely independent of inhomogeneous broadening and 
variations in matrix elements (Fig. |Tl]). 

Such multilevel molecular system can be described in terms of an effective two-by-two 
Hamiltonian Diagonalization of this Hamiltonian (Appendix y) allows us to analyze 
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FIG. 10: (a) Figure depicting resonant Raman experiments at 242 nm in which the 
molecule is excited to a high lying electronic state while fluorescence occurs at longer 
wavelengths, (b) By utilizing molecular coherence it is possible to enhance the rate of 
production of Stokes signal and mitigate fluorescence, (c) Scattering of the high frequency 
field at A3 off surfaces and other states can also be a source of confusion and this is 
eliminated by using a two photon IR drive. 
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FIG. 11: Energy level schematics for a three- level system to generate maximum coherence 
between the levels \b) and |c) via fields £i and £2- These fields are off-resonant with the 
electronic detuning 6lo and possibly also with the Raman detuning Auj which can vary 
in time, thus chirping the pulses. After preparing the coherence pbc with fields £1^2, the 
probe field £3 gives rise to the anti-Stokes field £4. 

the evolution of the system by drawing analogies to 2-state systems. If the excitation is 
applied resonantly (Au; = 0), such that the initial state of the system (the ground state 
|c) is projected onto the new basis formed by the eigenvectors |+) and |— ) [Eq. ( |C3|) ], the 
system undergoes a sinusoidal Rabi flopping between states \b) and |c). In this situation one 
can choose to apply a 7r/2 pulse in order to create the maximal coherence \pbc\ = 0.5. 

Alternatively, one can apply an excitation at a finite detuning Au, to allow all popula- 
tion, which is initially in the ground state, to follow the eigenstate |+) adiabatically. The 
coherence pbc is then 

1 



Pbc 



sin 9 e'^P 



(4) 



For molecular systems with large detunings, the Stark shifts A and D are approximately 
equal and 9 = tan^^ (2|i?|/2Aa;), where B is the effective Raman Rabi frequency (Appendix 
0). One method of achieving the condition \pf,c\ = 0.5 is to choose Au and to increase 
the product of the two incident fields until 9 is near 90°. This is done adiabatically with 
the product of the fields changing slowly as compared to the separation of the eigenvalues. 
Instead, at a fixed field, one may allow Au; to chirp from an initial value toward zero. 

We note that earlier, Grischkowsky |^ and Oreg et al. |^ have described preparation 
mechanisms in two-state and multi-state systems, and Kaplan et al. [H] have predicted 
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existence of 27r Raman solitons. 



2. Fractional STIRAP 



In an all-resonant A scheme (Fig. |T2|, with 6uj = Aoo = 0) maximal coherence can be 
prepared between the levels b and c in a fractional stimulated Raman adiabatic passage 



(STIRAP) set up by a counterintuitive pulse sequence 0, such that the population 

of the upper state a is always zero and fluorescence from this state is eliminated. This can 
be accomplished via a counterintuitive sequence of two pulses at frequencies Uab and Uac- 
Under the condition of adiabatic passage, the molecule in the initial state \b) is transformed 
into a coherent state {\b) — \c))/\^. 




FIG. 12: Energy level schematics for the generation of maximum coherence between the 
levels 1 6) and |c) via fractional STIRAP by counterintuitive pulses 1 and 2. After a time 
delay of r the pulse resonant with \a') — > \b) transition produces a signal at uja'c- 

The principle behind a STIRAP process is the adiabatic theorem as applied to the time- 
varying Hamiltonian H{t). If the system at time to is in an eigenstate of H{tQ), and the 
evolution from to to ti is sufficiently slow, then the system will evolve into the eigenstate of 
H(ti). The three-level atomic system driven by two fields has three eigenstates, one of which 
is a linear superposition of only the lower levels h and c. The time dependent amplitudes of 
this eigenstate depend on the pulse shapes of the fields at frequencies Uab and ujac- Thus, 
by an appropriate pulse shaping, it should be possible to prepare a maximahy coherent 
superposition of states h and c as shown in Fig. |l^. The expressions for the Hamiltonian 
and the corresponding eigenstates are given in Appendix 0. 

Comparing different schemes for the preparation of maximal coherence, we note that 
the required laser power is much lower for the all-resonant scheme, but in the case of 
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biomolecules, UV lasers are required. The far-detuned scheme will work with more powerful 
infrared lasers, up to the point of laser damage. As for the comparison of adiabatic and 
non-adiabatic regimes, we should note that the adiabatic scheme may turn out to be more 
robust, because it does not rely on a particular pulse area and works for inhomogeneous 
molecular ensembles and non-uniform laser beams. 



3. Femtosecond Pulse Sequences 



In a series of beautiful experiments K. Nelson and coworkers have generated coherent 
molecular vibration via a train of femtosecond pulses, see Fig. 0. They nicely describe their 
work as: "Timed sequences of femtosecond pulses have been used to repetitively "push" 
molecules in an organic crystal . . . , in a manner closely analogous to the way a child on a 
swing may be pushed repetitively to reach oscillatory motion." 

An interesting aspect of this approach is the fact that the individual pulses need not be 
strong. Only the collective effect of many weak pulses is required. This may be helpful if 
molecular "break-up", due to strong £i and £21 is a problem. This will be further discussed 
elsewhere. 




FIG. 13: Femtosecond pulse train prepares molecular motion, described by Qbc- The 
molecular state is then probed by field 3 so that field 4 is produced. 



B. Adaptive Evolutionary Algorithms 

So far we described how one-photon and two-photon resonant pulse sequences can be 
used to produce a coherent molecular superposition state. The idea is that once this state 
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is created, a delayed pulse can be applied in order to produce Raman scattering which will 
bear the signature of the molecular system. The Raman signal is expected to be optimized 
when the molecular coherence is maximal. 

In general, however, things are complicated by the Franck-Condon factors. As an example 
of a more complicated situation consider Fig. |T^. There we see a multilevel system with the 
ground state |c) and the next state \d). Matrix elements with Franck-Condon overlap factors 
yield a weak transition between \d) and \a) or \a') as indicated in the figure. However, given 
an appropriately nonlinear ground state potential it is quite possible that the next vibrational 
state could be off-set and have the appropriate position of the peaks of the wave function 
in order to maximize the Franck-Condon overlap. In this way, the coherence between the 



states 1 6) and |c) of Fig. |T4| could still serve as a strong generator of the anti-Stokes radiation 
on the I a) to \b) transitions. 




FIG. 14: Schematics to generate coherence between the levels \b) and |c). The shaded 
curves represent the probability distribution of the interatomic distance for different vi- 
brational and electronic states. 



Given such a configuration, it is not hard to see how to prepare the ground state co- 
herence. We could, for example, radiate the molecule with a chirped Raman pair so as to 
generate maximal coherence between \b) and \d) and then follow that by a Raman vr— pulse 
transferring the population from \d) to |c). In this way, maximal coherence between |c) and 
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\b) would be prepared. And as indicated in the figure, strong matrix elements would be 
expected between \b), |c) and the first excited vibrational state of the electronic potential. 

The preceding example shows that even in a simplified ideal few-level system, preparation 
of maximal molecular coherence may require application of a complicated pulse sequence. 
For large bio-molecules the level structure is not only much more complex, but usually 
unknown. We now consider how search algorithms can be used to find the optimal pulse 
sequence for a complicated molecule with an unknown Hamiltonian. This approach will 
eventually lead to an efficient generation of "molecular fingerprints" . 

In order to achieve this goal we will need to (1) utilize a technique for preparation of 
complex shaped pulse sequences; (2) find the particular pulse sequences, required for the 
excitation of the particular bio-molecules and for the production of spectral signatures, 
which will allow one to distinguish (with certainty) the target biological agent from any 
other species. 

Pulse shaping techniques already exist; they are based on "spectral modification". First, 



a large coherent bandwidth is produced by an ultra-short pulse generation technique [48 



Then, the spectrum is dispersed with a grating or a prism, and each frequency component is 
addressed individually by a spatial light modulator (a liquid crystal array |^ or an acous- 
tic modulator This way, individual spectral amplitudes and phases can be adjusted 
independently. Finally, the spectrum is recombined into a single beam by a second disper- 
sive element, and focused onto the target. This technique allows synthesis of arbitrarily 
shaped pulses right at the target point, and avoids problems associated with dispersion of 
intermediate optical elements and windows. 

A particular shaped pulse sequence can be represented by a three-dimensional surface 
in a space with frequency-amplitude-phase axes. Each pulse shape, which corresponds to a 
particular 3-D surface, produces a molecular response. The problem is to find the optimal 
shape. The search space is too large to be scanned completely. Besides, many local optima 
may exist in the problem. The solution is offered by "global search" algorithms (such 



as adaptive evolutionary algorithms) ||5T|, In this approach the experimental output is 
included in the optimization process. This way, the molecules subjected to control, are called 
upon to guide the search for an optimal pulse sequence within a learning loop |T3| . With the 
proper algorithm, automated cycling of this loop provides a means of finding optimal pulse 
shapes under constraints of the molecular Hamiltonian and the experimental conditions. No 
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prior knowledge of the molecular Hamiltonian and the potential energy surfaces is needed 
in this case. 



This adaptive technique was developed for coherent control of chemical reactions p2 
The idea is that the pulses can be optimized to produce desired chemical products. In our 
problem we want to optimize Raman generation. In this case both preparation and reading 
pulses can be adaptively shaped in order to maximize the signal. Fig. ^ shows schematics 
for the experimental setup that implements these ideas. 
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FIG. 15: (a) Experimental setup for the implementation of adaptive techniques; (b) 
Figure depicting amplitude of possible optimized Raman preparation pulse sequence 8182- 
Not indicated is the fact that the "reading" pulse £3 can also be profitably considered as 
a learning algorithm variable. 



Generated spectra will be different for different molecular species. And our task is not 
only to maximize Raman generation, but also to identify spectral patterns characteristic of 
particular species and maximize the difference in the spectrum produced by the target bio- 
molecule from spectra produced by any other bio-molecules. The key idea here is to apply 
the same adaptive algorithms in order to learn these optimal "molecular fingerprints" . 

We note that the complexity of the molecular level structure is not so much a problem as 
a solution to a problem. We can take advantage of the richness of the molecular structure, 
and the infinite variety of possible pulse shapes, in order to distinguish different species with 
the required certainty. 
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V. POSSIBLE FAST CARS MEASUREMENT STRATEGIES FOR DETECTION 

OF BACTERIAL SPORES 



Having presented the FAST CARS technique in some detail we now return to the question 
of its apphcation to "fingerprinting" of macromolecules and bacterial spores. Some aspects of 
the technique seem fairly simple to implement and would seem to hold relatively immediate 
promise. Others are more challenging but will probably be useful at least in some cases. 
Still other applications, e.g., the stand-off detection of bioaerosols in the atmosphere present 
many open questions and require careful study. In the following we discuss some simple 
FAST CARS experiments which are underway and/or being assembled in our laboratories. 




Hang off to fuiTher 
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FIG. 16: (a) Two-stage detection scheme. Airborne particles are irradiated with a UV 
laser. Those particles which fluoresce are selected for further testing, (b) FAST CARS 
testing of the preselected particles is like "singing a song" and listening for the standing 
ovation after an intricate aria. That is, the signal is generated only when a particular 
molecule hears its molecular melody. 



A. Preselection and hand-off scenarios 

At present, field devices are being engineered which will involve an optical preselection 
stage based on, e.g., fiuorescence tagging. If the fiuorescence measurement does not match 
the class of particles of interest then that particle is ignored. When many such particles 
are tested and a possible positive is recorded, the particle is subjected to special biological 
assay; see Fig. |T6|a. Such a two stage approach can substantially speed up the detection 
procedure. The relatively simple fiuorescence stage can very quickly sort out many uninter- 
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Conducting surface 
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FIG. 17: Surface Raman payoff: Proof of principle we are carrying out, utilizing a cold 
(He^ temperature) sample so as to enhance the Raman signatures and minimize 7;,c- 



esting scattering centers while the more sophisticated Raman scattering protocol will only 
be used for the captured "suspects" . 

The properly shaped preparation pulse sequence will be determined by, e.g., the adaptive 



learning algorithm approach as per section P^VB| . The amplitude and phase content of 
the pulse which produces maximum oscillation may be linked to a musical tune. Each 
spore will have a song which results in maximum Raman coherence. A correctly chosen 
"melody" induces a characteristic response of the molecular vibrations — a response which 
is as unique as possible for the bacterial spores to be detected. Playing a melody rather 
than a single tone is a generalization that enables us to see a multidimensional picture of 
the investigated object. We note that the optimization can (and frequently will) include not 
only the preparation pulses 1 and 2 (see Fig. |T^), but also the probe pulse 3, in particular, 
its central frequency and timing. Analysis of the response to such a complex input is a 
complicated signal processing problem. Various data mining strategies may be utilized in a 
way similar to speech analysis. 

However, taking into account the fact that we work with femtosecond pulses chained in 
picosecond to nanosecond pulse trains, the whole analysis can be very short. In particular, 
if we recall the long sampling time of the complete fluorescence spectra of being ^ 15min, 
our estimation of a microsecond analysis is a very strong argument for the chosen approach. 
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FIG. 18: Influence of various radicals on the frequency of the symmetric NO2 vibrations. 



B. Possible further Raman characterization 



After a suspect particle has been targeted, it may be subject to a whole variety of inves- 
tigative strategies. Raman scattering off a flying particle can be very fast, but not necessarily 
the most accurate method. It will be very useful to pin the particle on a fixed surface and 
cool it down to maximize the decoherence time T2 so that the characteristic lines are nar- 
rowed down. The particle can be deflected by optical means (laser tweezers, laser ionization, 
etc.) and attached to a cooled conducting surface (see Fig. |l^). Cooling to liquid helium 
temperature would enable us to enhance the dephasing time from T2 ~ 10~^^sec at room 
temperature to T2 ^ 10~^sec at a few degrees Kelvin. 



C. Possible spore specific FAST CARS detection schemes 

We conclude with some speculative observations for long range (stand-off) measurements. 

The chemical state of DPA in the spore is of special interest to us because the stuff we 
hang on the DPA molecule will determine its characteristic Raman frequency. To this end, 
we quote from article |5^ by Murrell on the chemical composition of spores: "When DPA 



is isolated from spores it is nearly always in the Ca-CDPA chelate but sometimes as the 
chelate of other divalent metals [e.g. Zn, Mn, Sr etc.] and perhaps as a DPA-Ca amino 
complex." 

Thus, since each different type of spore would have its own unique mixture of metals and 
amino acids, it may be the case that the finer details of the Raman spectra would contain 
spore specific "fingerprints." This conjecture is supported by Fig. where the difference 
between the DPA Raman spectra of the spores of Bacillus cereus and Bacillus megaterium 
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FIG. 19: Nuclear magnetic resonance ^^C Cross polarization magic angle spinning (CP- 
MAS) spectrum of (a) outer coatless mutant Bacillus subtilis 322 spores and (b) dormant 



Bacillus cereus spores. Complete spectra is found in |54| 



IS encouraging. 

The open question is: to what extent is the DPA Raman spectra sensitive to its environ- 
ment? That we might be able to achieve spore specific sensitivity is consistent with the well 
known fact that substituents, e.g., NO2 experience a substantial shift of their vibrational 



frequencies when bound in different molecular configurations, see Fig. Furthermore, 
recent NMR experiments show spore specific fingerprints due to the local environment 
(see Fig. 

Clearly there are many oportunities and open questions implicit in the FAST CARS 
"molecular melody" approach to real time spectroscopy. However it plays out, this com- 
bination of quantum coherence and coherent control promises to be a fascinating area of 
research. 
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APPENDIX A: THREE LEVEL SYSTEM DRIVEN BY TWO FIELDS 

In this appendix we consider the density matrix approach to the Raman scattering and 
discuss various hmiting cases. In particular we present a semiclassical treatment in which 
the field evolution is described by Maxwell's equations and the atomic system by the density 
operator. 

We consider a three level atomic system in the A configuration with upper level a and 
lower levels b and c. The a — > 6 transition is driven by a field at frequency i>i and the a — c 
transition is coupled via a signal field at frequency z/2. 

The Maxwell's equations lead to the following equation for the evolution of the signal 
field, 

- -^T^] E = l^oP, (Al) 



where E is the electric field vector of the emitted light, /lo is the vacuum permeability, and 
P is the medium polarization. We write the electric field using the slowly varying envelope 
E as 

E = ^g{z,t)e-'^'''-'''+'^^ + C.C. (A2) 

and write the polarization in terms of the slowly varying quantity P as 

P = lp(2,t)e-^(^*-'^^+'^) + C.C. . (A3) 

Working within the slowly varying amplitude and phase approximation, Eqs. ( [AID , ( |A-2[ ), 
and (|A|) yield 

Ids dS 1 ^ ,^ 

where 

V = 2Af{z,t)p,,gabe'^^'-'^~"^\ (A5) 

Here pab is the dipole moment matrix element, Qab is the off-diagonal element of the density 
matrix for the molecular levels a and b, and N'iz, t) is the volume density of the molecules. 
We note that the electromagnetic field is determined by Qab- 
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FIG. 20: Schematic diagram of the Raman interaction for the three-state system. The 
driving fields 1 and 2 are in Raman resonance with the two lower levels b and c, but 
are generally off-resonant with the transitions to the excited state a (with detuning A). 
Decoherence and decay of the optical transitions ab and ac are characterized by the rates 
Jab and 7ac, the spontaneous decay rate of level 6 to c is 71, and the decoherence rate of 
the be transition is j^c- 



Here we recall the main results of |40] and apply them to our physical situation (see Fig 
20]) . First we define the decay rates 

7 + 7i 



lab 

Icb 



abi 



1 

2 

71 



+ 7^ 
2 ' leal 



+ llb- 



(A6) 
(AT) 
(A8) 



Here the decay rates from a to 6 is 7^, and from a to c is 7c, and 7 = 7fe+7c. Population decay 
rate from 6 to c is 71. Purely phase decays are designated by the superscript p. Complex 
dephasings are defined as 

Pafe = lab - i'^ab, (A9) 
Pea = lea + iKc, (AlO) 
Tcb = Icb, (All) 

where Aac = uJac — ^1 and Aab = 1^2 — ^ab- In the following we assume Raman resonance, i.e.. 



are 



-Aab = A. The main working equations for the off-diagonal density matrix elements 



Qab = -TabQab + i^2{Qbb " Qaa) + i^lQcb, 
Qca — ^ caQca ~\~ ^^l(f?aa Qcc) ^^2&ti; 
Qcb = -TcbQcb - i^2Qca + i^lQab- 



(A12) 
(A13) 
(A14) 
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The equations for the populations are 



Qbb = i^lQab - i^2Qba + IbQaa - llQbb, 
Qcc = i^lQac - i^lQca + IcQaa + llQbb, 



(A15) 

(A16) 



and Qaa is obtained from 



(A17) 



Here Qi^2 are the Rabi frequencies of the fields having frequencies z/i 2, i.e., 

Pac^l 



^2 



h ' 

Pab£2 



(A18) 
(A19) 



The steady state solution of these equations can be obtained by setting all time derivatives 
equal to zero. The result is 



Qab = -T^[^2{Qaa - Qbb) " ^iQcb], 
i ab 



(A20) 



and 



Qca = [{Qaa " Qcc) ab'^ cb + \^l\^) + (Qbb " ^'aa)|J^2p] 



Qcb 



V 



[{Qaa - Qbb)^ca + {Qaa " Qcc)^ 



ab\ 



(A21) 
(A22) 



where the common denominator is 



ca^cb + '^ab\^2\ +rca|fil| 



(A23) 



Equation ( |A2U| ) is our main working equation. The first term in the parenthesis of Eq. 
( |A2(J| ) is responsible for the stimulated emission or absorption on the \a) ^ \b) transition. 
The second, Qcb, term describes the Raman conversion of z/i into z/2. As we can see, it is 
crucial to have the element Qcb as big as possible for optimum Raman conversion. In the 
following we present the main properties of Qcb for various experimental and conceptual 
configurations and consider various limiting cases for the Raman processes. 
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1. OfF-resonant Raman process, weak driving 

First we consider the Raman effect far from electronic resonance (A ^ 7a6,7ac)- In this 
case the upper level a is almost completely depopulated and can be eliminated from the 
dynamics. The system behaves basically as a two-level system with the Rabi frequency of 
oscillations between levels h and c equal to Vt^f = |^]*^72|/A. For dephasing rate ^ ^r, 
we find the steady state solution for the lower state coherence Qcb as 

lim ^cfe = i^-^iQcc - Qbb), (A24) 

which, after substituting into Eq. ( |A2CI| ) , yields the density matrix element responsible for 
the U2 radiation as 



Raman ■ ^2 



Qab 



lab - A 



|2 



[Qaa - Qbb) - ^^^^{Qcc - Qbb) 



(A25) 



7fecA 

Note that for the amplification of the z/2 field the so-called Raman inversion g^c — Qbb is 
necessary. In the simplest case when almost all the molecules are in the lowest state c, 
~ 1, Qbb^O, Qaa^O, Eq. (|A25|) reduces to 



= (A26) 

A A'jbc 

2. Off-resonant Raman process, strong driving 

In the case of strong driving the transition between states b and c can be saturated so 
that both these levels are significantly populated (with virtually no population in the excited 
state a) . In particular, one obtains from Eq. ( |A15| ) the steady state populations of the lower 
levels 

in^pab - i^2pba I . „„s 

pbb — , l^^'J 

7i 



and Pec = 1 — Pbb- Substituting Eq. ( |A20|) into Eq. (|A27] ), and assuming, without loss of 
generality, that the Rabi frequencies for the laser fields are real {Vti = VL\,VL2 = we find 
that for the far off-resonant case {Tab ~ ^A) 

Pbb = (A28) 
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Rearranging Eq. ( |A22|) and inserting Eq. ( |A28| ), under the condition 



Ibc > 



+ 



we find that 



Pbc 



iA iA 



iATbc 'yiA^Tbc 
Solving Eq. (|A30|) for pbc we obtain 



(A29) 



(A30) 



Pbc 



76c + 



A V'" ' 7i^^ 

It can be shown that the maximum magnitude of the coherence \pbc\ is given by 



\Pbc 



71 
I67 



be 



This maximum occurs when 



^1^2 _ a/tF)^ 

A ~ 2 ■ 



(A31) 



(A32) 



(A33) 



3. Resonant Raman process, weak driving 

In the resonant case (A = 0), the Raman effect has much in common with the scheme of 
lasing without inversion (LWI; see, e.g., The lower state coherence is then found to be 

[{Qaa - Qbbhca + {.Qaa - Qc,^la}}\ 



Qcb ^ ^72^1" 



lablaclcb + lca\^l? + lab\^2\'^ 

Polarization responsible for the V2 radiation is then governed by 



.^2 

I — 

lab 



yQaa — Qbb) 1 — 



\^l\'^lab 



where D = -fablcajcb + 7afe|^2p + 7ca|^^lP• 
In the weak driving limit 2 ^ jab, lac, Ibc), almost all the population remains in state 
c (f?cc ~ 1, Qbb ^ 0, f?aa ~ 0). The atomic coherence of the lower levels is 

^1^2 



(A34) 



(A35) 



Pcb -- 

ladbc 

and the corresponding expression for the density matrix element pab is 

.^2 W 



Resonant Raman 
Sab 



—l- 



lab laclbc 



(A36) 



(A37) 
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4. Resonant Raman process, strong driving 

In the resonant case we can find the steady state solution of a strongly driven three-level 
system (|^^i,2| ^ 7i, 76c) in the form of the so-called dark state. It is useful to consider two 
superpositions of the lower states \b) and |c) defined as 

iB> = ,A38) 

where = a/I^^iP + |^^2p- The state \D) is completely decoupled from the upper state \a) 
{H\D) = 0), and is called the dark state. If the molecule starts in any state different from 
\D), the fields Qi and ^2 will promote it to the upper state \a), which decays to the lower 
states by spontaneous emission. If the molecule is in \D) it stays there unchanged. In this 
way the dark component of the state increases and finally the system will end up completely 
in the state \D). The populations of the levels a and b are paa = and pbb = l^iP/^^^, and 
the dark state coherence of the lower states is, as given by Eq. ( |A39| ), 

^cb = (A40) 

Maximum Raman coherence is achieved with = \Q2\ so that Ip^'^^l = 1/2. The element 
Qab [see Eq. ( |A2(]| )] responsible for the radiation is given by 



Dark Raman 

lab 



^2{Qaa - Qbb) + 



0, (A41) 



9? 

i.e., there is no radiation from the dark state. To obtain Raman signal one can, e.g, switch 
off the field ^2 after reaching the maximum coherence \q^^\ = 1/2. Then the density matrix 
element Qab goes as 

Max coherence ^ ^1 ^ A 

27a6 

To summarize the main results of this appendix, we display the density matrix elements 
responsible for the signal field generation in Table |. 
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APPENDIX B: COOPERATIVE SPONTANEOUS EMISSION 

Since the DPA is contained in the small {V ~ (1/^)^) volume of the core and since the 
number of participating DPA molecules is huge (A^ ^ 10^) the possibility of Dicke superradi- 
ance should be carefully considered. At the outset, we should note that the short dephasing 
times of vibronic levels (T2 is in the picosecond range) and the unknown inhomogeneous 
dephasing time T2* will tend to wash out the effect. Still, if the spontaneous emission time 
(r ~ 1 — 10 nanoseconds) is reduced by even a small fraction of N, then cooperative spon- 
taneous emission may be important. 

Motivated by the proceeding we now turn to a short review of the effect from the present, 
FAST CARS, perspective. We focus on two levels \a) and \b). The previous pulses ei and 
£2 beat the coherence between \b) and |c) so that at time t = 

\^m = B\b) + C\c), (Bl) 

the third pulse, arriving at time T, promotes |c) — > |a) so that at time t = t 

{"^{t)) =A\a)+B\b) + C'\c). (B2) 

Hence the atoms at time r are in a coherent superposition \a) and \b) and this affects the 
rate of spontaneous emission from the N molecules involved. 

For the present purposes it is best to start with the Hamiltonian describing the interaction 
of N two level systems (|a) and \b)) and the quantized radiation field given by 

Hint = ^i,khga+aie'''-'' + adj (B3) 

where, as discussed, e.g., in p, (7^ is the coupling frequency between the molecules and the 
j^th jj^q(Iq Qf ^j^g field, is the usual creation operator and cxj = (|6)(a|)j is the lowering 
operator for the i*'' molecule. 

Consider next the case in which all molecules are located at = tq in a volume small 
compared to the radiation wavelength A, and noting that the coupling constant is a slowly 
varying function of k we may write 



9ke 



iko-ro 



37 



The important point being that the lowering operator described by Eq. ( |B4D is symmetric 
in the molecular lowering operator cTj. 

To glean the physics from Eq. ( |B4| ) it is enough to consider 3 molecules which are all 
started in their upper state, so that 



\a1a2a3). 



(B5) 



The state (p5| ) evolves under the influence of the fleld according to the interaction ( |B4D , 
and the molecular state develops into |\l/3) al'^s) +P\'^2) where a and P are uninteresting 
constants and 

1^2) = J^o-jl^s) = \b1a2as) + \a1b2a3) + \a1a2b3). (B6) 

i 

In general the symmetric interaction (p^) only couples the 1010203) state to the symmetric 



states of Fig. 21 
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FIG. 21: Showing that the symmetric combinations of the molecular states |aj) and |6, 
1=1,2,3 are just the angular momentum states for j — 



_ ™ . 1_3 _|_ 1 

2 ' '"j ~ ^ 2 ' 2 • 



Thus we may, following Dicke describe our N "spins" by an effective angular momentum 
of magnitude j = |(iV+ + iV_) and projection nij = |(iV+ — A^_). 

Using this "angular momentum" picture and recalling that the matrix element governing 
the transition rate between state \j,mj) and \j,mj + 1) is given by 



^ = (j, rrij + l\J-\j, rrij) = J (j - mj)(j - mj + 1) 



(B7) 



we have two interesting limits. 

First the case when all spins are up i.e., = N, = (which corresponds to the case 
of all molecules in | a) ) ; then j — ^ 



N++N ^ 1 



2 and rrij = ^ — > y and we have 
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/i = (j, m,- + 1 1 jL I j, TJij) = Vn. (B8) 

Therefore the radiation rate goes as /i^ = A^, which is the normal result for N independent 
radiators. 

In the other case suppose that rrij = so that = y and = y. Then the relevant 
matrix element is 

/i = (j, -l\J-\j, 0) = V7(JTT) - N, (B9) 

and the radiation rate goes as N^. That is the state \j, 0) decays N times faster than the 
state and therefore the state for which nij = (equivalently = y and = y) 

is said to be "superradiant." In fact, any such state for which \mj\ <^ \j\ is said to be 
superradiant. 

Returning now to our problem in which the levels \a) and \b) are coherently prepared, 
suppose that A and B of Eq. (pl| ) are both An ensemble of N such coherently prepared 
molecules will be described by a superposition of rnj levels sharply peaked about rrij = 
and is thus superradiant. 

The major potential payoff from our perspective is the hope that the coherent resonant 
Raman process may direct the molecules from c ^ a — ^> 6 on a time scale which is of the 
order r/N where r is the spontaneous fluorescence time. This would be interesting since, 
as emphasized by Nelson and coworkers "This [utilization of resonant Raman spectra] 
is possible, however, only in the absence of fluorescence interference." 

But, as is discussed in Section there are many open questions. For example, will the 
coherence decay rates T2 and Tg* spoil the effect? These and other questions need to be 
addressed by careful experiments and further analysis. 
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APPENDIX C: EIGENSTATES AND EIGENVECTORS FOR A MOLECULAR 
SYSTEM WITH LARGE ONE-PHOTON DETUNINGS. 

In the case of large one-photon detunings 5uj (such that only small virtual excitation of 
the upper state is possible), the molecular system of Fig. ^ is described by an effective 
two-by-two Hamiltonian 



H - ^ 



A Q 
Q* D-2Auj 



(CI) 

where Q = VLi,VLf./5uj is the effective Raman Rabi frequency, A and D are dynamic Stark shifts 



see 



Il| for more details), and AVt is the (two-photon) Raman detuning. The eigenvalues 



of the Hamiltonian ( |C1D are 



-{A + D- 2Au) ±^{A-D + 2Auy + 4\Q\' 
E± = ^ , (C2) 



and the corresponding eigenvectors are 



In the case when A ^ D we have 



A- AujT ^/AuJ^^T~\^ 



2 

The probabilities of occupation of states \b) and |c) are 



(C4) 



|2 



+ (A + 2E±)2 







2 






2 + Acj2±Acuv/Acj2_^ 




1 



(C5) 



Pc ^ ' (C6) 



2 + Acu2 ± Auy/Auj^ + \Q\ 



Thus, if we start at t — >• — oo on eigenstate |— ) with Au — lAc^l, we have Pb(— oo) = 
and Pc(— oo) = 1. Chirping Au such that for t — > +oo we have Au \Auj\, the |— ) 
eigenstate changes adiabatically such that the probabilities get exchanged to Pfe(oo) = 1 
and Pc{oo) = 0. At the moment when Au = we have Pb = Pc = 1/2 and pbc = — l/2e*''^, 
where Q = Ifile**^. Thus, if we turn off the pulses at this time, we reach full coherence. 
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APPENDIX D: FRACTIONAL STIRAP 



We consider a three-level atomic system driven by resonant pulses at frequencies Uab and 



ujac with Rabi frequencies fii(t) = pab£-i{t)/h and Vt2{t) = pac^2{t)/Ji, respectively (se ^5 
for more details). The Hamiltonian for the system in the slowly varying amplitude and 
phase approximation is 

H = |(fii(t)|a)(6| + n2it)\a){c\ + H.c), (Dl) 

where H.c. means Hermitian conjugate and \a) is the excited state. This system has a 
dark state. The eigenvalue of the Hamiltonian for this state is equal to zero, = 0, i.e., 
H\D) = Xd\D) = 0. For the scheme shown, the dark state is 

^ Q^{t)\c) - Q,{t)\b) 

y/\n,itw + \n,{tW 

This state mixes the ground states \b) and |c) and is independent of the excited state \a). 

The sequence of pulses with Rabi frequencies fli{t) and ^l2{t) is such that, for the molec- 
ular system initially in the state \b), Qi{t)/Q2{t) ^Oatt—>- — oo and Qi{t)/Q2{t) tan^ as 
t — > +00. Thus the state evolves adiabatically from \b) to the coherent state cos0|6) — sin^|c). 
For the case when fii(t)/i72(^) 1) d = 7r/4 and the resulting state is — \c))/^/2. We 
can see that, unlike the case of conventional STIRAP, the two pump pulses vanish simulta- 
neously. In the conventional STIRAP the main aim is to bring a system from one low-lying 
energy eigenstate to another without populating an excited state on the way. The pump 
pulses in the conventional STIRAP are timed such that the pulse connected to the target 
state starts before the beginning and ends before the end of the pulse connected to the initial 
state (i.e., the so called counterintuitive sequence). The main aim of the fractional STIRAP 
is to bring the system from an energy eigenstate to a preselected superposition of two lower 
lying eigenstates, without populating any other state. This task requires a modified timing 
of the pump pulses. In contrast to the scheme discussed in the preceding Appendix |C|, the 
pump pulses of the fractional STIRAP are resonant which means that much weaker field or 
shorter duration of the pulses is required. 
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